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We analyze the effect of a form factor in the magnetic contribution to K'^ ■k'^-k'^^. We 
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I. INTRODUCTION 

Non-leptonic kaon decays are an iniportant tool 
to study weak interactions P, 0, H, H, H, @, 01 ■ Ra- 
diative non-leptonic kaon decays, such as K^^s ^ 
7r+7r~7 and —y 7r+7r''7 are dominated by long 
distance contributions. The study of these decays 
leads to chiral tests and in principle to disentangle 
the small short distance contribution. This small 
short distance part may lead also to interesting CP 
violating observables in the Standard Model (SM) 
and Beyond (BSM) 0]. 

K — > 7r7r7 amplitudes contain two types of con- 
tributions: inner bremsstrahlung (IB) and direct 
emission (DE). Due to the pole in the photon energy 
the IB amplitude is generally enhanced compared to 
DE; however the IB components of Kl 7r+7r~7 
and tt^tt'^j are suppressed due, respectively, 

to CP invariance and to the A/ = 1/2 rule. Then, 
the DE amplitude of these channels, that are the 
non-trivial part of these decays, might also be eas- 
ier to detect. DE contributions can be decomposed 
into electric (E1,E2) and magnetic (M1,M2) ones 

0. The magnetic contribution to 'k'^it~^ is 
accurately measured [l3, [Hf and a clear and large 
photon energy dependence has been found. 

The question of the presence of the form factor, 

1. e. the size of vector meson dominance in weak am- 
plitudes, is motivated not only by phenomenological 
reasons, as in Kl 7r+7r^7, but it has also theo- 
retical motivations For instance, it improves the 
matching of long and short distance contributions 
0. 

Published data for K^ n^n^j are consistent 
with a dominant magnetic amplitude and no evi- 
dence for El transitions [ll|. However, preliminary 
data from NA48/2 at CERN show a non- vanishing 
interference, due to El transitions [H, |T3|: the size 
of this contribution will shed light on 0{p'^) chi- 



ral perturbation theory (xPT) counterterm coefR- 
cients (T], [H, [l^. The energy dependence of the 
K^ —> 7r+7r'^7 magnetic amplitude has not been 
tested/observed yet; we think it is important to im- 
derstand if the form factor in the magnetic contri- 
bution affects the determination of the electric con- 
tributions: in this paper we complement Ref.[l6| 
with this perspective. 

Long ago N. Christ used a particular set of K^ 
TT^TT^j Dalitz variables: the photon energy and the 
charged pion kinetic energy in the kaon rest frame 
[13. NA48/2 [H, nil wants to perform an anal- 
ysis using these variables too; in the following we 
give a kinematical distribution in these variables ac- 
counting for the — rn„o corrections. Further- 
more, the results in Ref.[lal considered essentially 
the central value of branching ratio from E787 [131 , 
(4.7±0.8)xl0-6. Recent results [H [Hi seem to 
prefer smaller values, (2 ^ 3) x 10~^; while this 
is not shocking by itself, it has some impact if we 
include the form factor as we shall see. We dis- 
cuss K — > 7r7r7 kinematics and Low theorem in Sec- 
tion 2; in Section 3 we carry out some theoretical- 
phenomenological considerations on A' — > ■ktt^. Nu- 
merical results are summarized in Section 4. 



II. KINEMATICS AND LOW THEOREM 



The general invariant amplitude of K t:t:^ can 
be defined as follows [1, [l^j 

A[Kip) ^ 7ri(pi)7r2(p2)7(9,e)] = e^((z)A-f^(g,pi,p2), 

where e^{q) is the photon polarization and Mf^ is 
decomposed into an electric E and a magnetic M 
amplitudes as 
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The double differential rate for the unpolarized pho- 
ton is 
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where = m^^/mif. Low theorem establishes a 
precise relation between radiative and non-radiative 
amplitudes in the limit oi [2l|. Then 

we can generally write the relation between the 
bremsstrahlung amplitude, -E/s, in K — > -k-k^ de- 
cays and the on-shell amplitude, A{K tttt): 
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where Qi is the tt^ charge. Direct emission ampli- 
tudes arc defined by subtracting this contribution 
from the total amplitude. In the case 



Eib{K+) 



( 3e^2 



\2mKZ+Z2, 



(4) 



where A{K^ —> tt+tt ) ~ ^A2e and z+ refers to 
the charged pion. Using the experimental value for 
B{K^ TT+TT^) we obtain the branching ratios for 
the Inner Bremsstrahlung shown in TAB. HI 
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(2.61) X 10"* 


[0,80] 


(1.84)xl0"* 



TABLE I: Inner Bremsstrahlung 

In Kl 7r+7r"7, the most common variables 
are: i) the photon energy in the kaon rest frame 
E* , and ii ) the angle 9 between the photon and 7r+ 
momenta in the di-pion rest frame. The relations 



between E*, 



and the z, are: 



pion kinetic energy in the i^"*" rest frame T*, and 
u) = {q-pK){q-p+)/{ml+m\) ^. These two 
variables are related to the Zi by 
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The advantage of these variables is that, under the 
assumption of constant Ede and Mde, one can 
easily disentangle the different contributions of the 
IB, DE amplitudes, and interference term between 
IB and DE 
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where 



A = A{K+ tt+tt") 



(9) 



A. Kinematics in Christ's variables T* and Ej 

Motivated by the upcoming NA48 measurements 
we have also studied the Dalitz plot distributions in 
the kinematical variables T* and the photon energy 
in the kaon rest frame. El. In this way we extend 
the work by N. Christ [T3I by including the terms 
proportional to — m^o; in these variables the 
double differential rate for 
as 



TT+TT 7 is written 
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where [3 = ^l- Aml/{mj^ - 2m rE*). Then the 
differential rate is 
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Since three photons will be detected in the A'+ 
tt+tt'^j measurement, it is very useful to study the 
differential rate as a function of: i) the charged 



where Eq is the Tr'^-energy, 
r 2 2 



and the Inner Bremsstrahlung differential rate is 
written in terms of angle, 6, between the 7r+ and 7 
momenta in the kaon rest frame 
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not a CP eigenstate, Ml and El transitions arc CP- 
allowed in ^ while for K^^ CP-symmetry allows 
Ml transitions and inhibits El transitions. The 
0{p^) contributions to electric transitions can 
be parametrized as 
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where = 9 • 10 GcV is obtained from the 
A/ =1/2 contribution to Ks tttt at ©(p^) and 
N^j^'' is the relevant C'(p^) weak counterterm com- 
bination H, [n, m 



FIG. 1: Dalitz plot of Christ's variables T* and E* 
and contour plot of the Inner Brcmsstrahlung ampli- 
tude: the right corner (in red) has a greater density. 
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It is possible to pass from eq. (fTU)l to cq. 
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A contour plot of the IB amplitude in the Dalitz 
plot of Christ's variables T* and E* is shown in 
FIG.l. After integrating ea. pO)) on E*, one obtains 

{E*=T:+m^^) 
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The second line is the theoreticalprediction based 
on the factorization model (FM) [E H [H, S [11] 
parametrized by the a coefficient fc/. From nai've 
dimensional analysis (NDA) we expect n'^^^ of order 
one; in fact this value is expected from VMD and 
factorization [ll,[il]. The sign in and HH) leads 
to a constructive interference among Ei and Ejb 
[2[i | but so far there is no experimental evidence of 
such interference, up to the new NA48 result [H, 
as we shall see. 

Actually, the present bounds [isl. Il9l| a re very 
close to the theoretical predictions (ll. Il4l. [T5| : in 
TAB. [n] we report instead the results for a de- 
structive interference, N^^'^ = 0.4, as reported by 
NA48 [13. The interference term scales linearly 
with n]^'^, then branchings for different n''^'^ val- 
ues are easily obtained. The same counterterm co- 
efficient, N^^'^, appears in the DE component of 
Ks 7r+7r~7 |20|. [23 |. however the present exper- 
imental bound [2J| is not at the level to compete 
with the one from 7r+7r"7. 

No El (CP violating) transitions for Kl have 
been observed yet [To| ; while the magnetic contribu- 
tions, Ml is responsible for the observed B(Xi 
TT+TT~j)]jE- The leading order magnetic amplitude, 
M{zi), to eq. ID), starts at 0(p*) 
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III. EXPERIMENTAL STATUS AND 
THEORETICAL PREDICTIONS 

Electric and magnetic dipole (El and Ml) transi- 
tions appear already at 0{p'^) [i3j[3- Since AT"*" is 



The subscripts L and + denote n^n^^ 
and tt'^tt^j respectively. The a^'s parts 

of the above amplitudes come from the local weak 
lagrangian Cf^^^ and are expected of order one 
[j^[23; these are also called direct contributions. 
They were originally derived using the factorized 
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form, current-current, of the weak chiral lagrangian 
with one current 0{p'^) from the WZW lagrangian 
and the other current from the usual 0{p) current 
[2^. Later it was found that also vectors and ax- 
ials (VMD) were contributing to these coefficients 
(oi's) [Hi. The factor "2" in ^ is the contribu- 
tion from the A'+ — ^ "7r+"^ 7r+7r°7, where we 
have first a pure weak transition and then a Wess- 
Zumino-Witten one and thus it is completely pre- 
dicted (indirect contribution). There are too many 
Oj's to be fixed phenomenologically and possibly 
large 0{p^) corrections to Kl 7r+7r~7 in (fT5|) 
from ry'-exchange. 

An interesting progress to the understanding of 
these decays has been driven by KTeV; in order to 
fit the Kl 7r"'"7r~7 data showing a clear photon 
energy dependence in the magnetic term, M{zi), 
KTeV has used a linear fit, a quadratic fit or a pole 
fit of the form 
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Interestingly, data prefer this pole fit to linear and 
quadratic fit [lol |. The rate and the photon en- 
ergy spectrum fix IHmI = (9.4 ± 0.8) • 10^'' and 
b = -1.243 ± 0.033 This phcnomenological 

description has been comforted by NA48 [IHl and 
KTeV ^ in the channel Kl Tr+Tr-e+e". Com- 
parison of ([T5)) with leads to the following the- 
oretical consequences: i) the value of hm fixes the 
size of the a^'s: ^ C(l) and ii) the presence of 
a relatively large component, b, with form factor 
tells us that VMD plays a major role in these co- 
efficients. The usual vector formulation, "y'"', is 
very successful to this description: there are vector 
and axial contributions to the a^'s [l5| . However 
the tensor formulation of the vectors, "]/^''", very 
successful in the strong sector, does not have any 
vector contributions to the ai's The lack of 

vector contributions to the Oj's and possible large 
0{p^) contribution to (fT3|) from jy'-exchange docs 
not explain the observed large form factor in (|15p 
and puts tension in the chiral expansion with large 
0{p^) and smaU ©(p") contributions Q. Actually, 
in Ref. [2^ it was already shown that while the 
odd-parity couplings to y — *■ P7 decays, relevant 
to the anomalous K —^ ttttj decays, had the proper 
QCD behaviour if we use the usual vector formu- 
lation, "y*", this was not the case if we use the 
tensor formulation of the vectors, "X/m""^ As a re- 
sult we believe in the large VMD contribution to 
the Oj's. 

We show in TAB. HI] an updated K+ tt+tt^j 
experimental and theoretical situation. We write in 
the first two lines the PDG 06 value along with its 



REPS. 


B{K+ ^ 7r+7r°7)DE 
T* e [55, 90] 


BNL E787 [18] 


(4.7±0.9)xl0-'' 


PDG 06 [llj 


(4.4±0.7)xl0"^ 


KEK-E470 [19] 
NA48/2 [131 


(3.8±0.8±0.7)xl0-'* 
(2.22 ± 0.13 ± 0.05) X 10"^ 


NA48/2 analysis [IS] 


T; e [0, 80] MeV 


B{K+ ^ 7r+7r«7)™T 


(-4.91 ±2.00) X 10"" 


B{K+ ^ 7r+7r«7)°'= 


(6.16 ± 0.79) xlO"'* 




Theory 
-range in MeV 


predictions 

B{K+ ^ ^+7r°7)iNT 


[55, 90] 
[0, 80] 

T^* -range in MeV 


-(3.52)xl0~'^ 
-(4.70) X 10"" 

B{K+ ^ 7v+n°-f);;;° 


[55, 90] 
[0, 80] 


(3.55) X 10"^ 
(6.57) X 10"'' 



TABLE IP Experimental and theoretical status 
The table shows in the first two lines the PDG 06 
value along with its most precise measurement, BNL 
E787 and in the next two lines two subsequent mea- 
surements, KEK-E470 and NA48/2, pointing towards 
a smaller value of the branching ratio. All these val- 
ues are obtained with vanishing interference. To com- 
pare NA48/2 with other experiments we extrapolated 
the NA48/2 value, obtained in the kinematical range 
< r* < 80 MeV, to the kinematical range 55 < T* < 
90 MeV, assuming a constant magnetic term. We report 
in the fifth and sixth row, the interference and the di- 
rect emission contributions determined simultaneously 
by NA48/2 The INT theoretical branching ratio 

with a value of the weak counterterm combination in 
eq. so to match the NA48/2 resuh, N'-^^ = 0.4, is 
then shown. Por comparison the DE theoretical branch- 
ing ratios obtained with ai = are reported in the last 
rows. 



most precise measurement, BNL E787; in the next 
two lines two subsequent measurements, KEK-E470 
and NA48/2, that as we can see, point towards a 
smaller value of the branching ratio. All these val- 
ues are obtained with vanishing interference. Also 
to compare NA48/2 with other experiments we ex- 
trapolated the NA48/2 value, obtained in the kine- 
matical range < T* < 80 MeV, to the kinemati- 
cal range 55 < T* < 90 MeV assuming a constant 
magnetic term. Interestingly NA48/2 has also done 
the analysis to determine simultaneously both the 
interference and the direct emission contributions 
[ill; we show in TAB.HJthis NA48/2 analysis show- 
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ing non-vanishing values for both the interference 
and the direct emission contributions. 

We also show in TAB. for comparison some 
theoretical predictions for the interference term and 
the direct emission term. Regarding the interfer- 
ence term we use a specific value of the weak coun- 
terterm combination in eq. (Ha): N'j^j = 0.4, so 
to match the NA48/2 result. Regarding the DE 
contribution we take as comparison the = 0. 

Following Ref. a correlated analysis of the 
— > 7r+7r''7 and — > 7r"'"7r~7 decays can be 
performed. VMD and phenomenology, i.e. eq. (fT5)). 
imply the following decomposition: 



M+(z,) = i/p°'" + A/: 



(16) 



It is interesting that all the weak couplings involv- 
ing vectors, and consequently M_|_ and M^^^ , may 
be written in terms of only one coupling, r/v [la]- 
Thus we parametrize all our ignorance in M^^^^ with 
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FIG. 2: Values of and rjv on the central solid line 
generate the E787 experimental value of the branch- 
ing ratio in TAB. |II] while the other lines are the 
borders of strips corresponding to la, 2a and 3cj 
deviation from the E787 central value in TAB. |ll] 
As a result the lower solid line corresponds to the 
central NA48 value in TAB. HH For r^v = and 
A+ = 2, the amplitude is dominated by the WZW 
pole. 
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is determined by the VMD couplings [1, [TH| fv ~ 
0.2 and hy ~ 0.037, well known phenomenologi- 
cally. The rest is written as a constant contribu- 
tion. 



where the parameter 
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In fact, recent data, shown in TAB. [Til point to- 
wards smaller values of the branching. 

NA48, in their preliminary analysis [l^, are able 
to study this decay in a new kinematical region: 
the charged pion kinetic energy ranges from to 
80 MeV; as we can see from FIG. 1, this region is 
more sensitive to the DE component. Furthermore 
assuming a constant M{zi) in ([T]) they find a non- 
vanishing i51-contribution. 

This result is very interesting; in particular 
though the size is comparable to the theoretical ex- 
pectations in ^ and TAB. IT] [J, [13, [H , the sign 
is opposite as discussed before. In the next section 
(19) we investigate the consequences of the presence of 
the form factor on this NA48 measurement. 



~ 2-|-3(2a3—a2)non-VMD+ other contributions. 

must be determined phenomenologically. (2a3 — 
fl2)non-VMD IS a non-VMD contribution to M~^; 
thus the experimental B{K'^ 7r"'"7r'^7)£)^ can be 
obtained by varying the two unknown constants 
(A+,77y). In FIG. 2, we vary respectively la, 2a 
and 3a deviations around the BNL E787 result [Tst 
similar to the PDG rcsidt in TAB. |TT1 So we can 
account properly the sensitivity of the W— and 
T*— spectrum from these phenomenological param- 
eters which we will discuss in the next paragraph. 



IV. EFFECT OF THE FORM FACTOR ON 
THE DIRECT EMISSION AMPLITUDE 



The target of the following numerical studies is to 
understand how to distinguish a constant magnetic 
term, practically rjv = 0, from the one with the 
form factor in eqs. JTC]) and (H?]), in FIG. 2. A 
related question is, as we shall see, if the possible 
presence of a form factor affects the determination 
of the interference term or even if the presence of the 
form factor could mimic the experimental evidence 
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of the interference term. It is important, we think, 
to quantify this effect. 

Previous experiments were able to study only the 
kinematical range T; S [55, 90] MeV recen- 
tely NA48/2 has been able to uncover almost the 
full range T* € [0,80] MeV [H. The region with 
small values for T* is more sensitive to DE tran- 
sitions. Indeed the explicit IB Dalitz plot contour 
in the T* and W variables generates a plot similar 
to FIG. 1. Due to the sensitivity of the kinemat- 
ical distributions from the size of DE branchings 
we will discuss both possibilities Bexp(E787) and 
Bexp(NA48) in both kinematical regions. 
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the central (lower) full line of FIG. 2 for the upper 
(lower) plot. 

Then we plot the W^-spectra with T* G [55,90] 
MeV in FIG. 4 and T* G [0, 80] MeV in FIG. 5. In 
each case we consider form factors with rjv = 0.5 
(dashed curves) and rjv = 1-5 (dotted curves); the 
corresponding values of are determined from 
FIG. 2. 
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FIG. 3: Normalized T*-spectra (T* G [0, 80] MeV) 
for the DE magnetic amplitude with the E787 
branching ratio [ia ] in TAB. |lT] (upper plot) and 
the NA48/2 branching ratio in TAB. [6] (lower plot). 
The solid curves correspond to a constant ampli- 
tude, the dotted curves to a magnetic form factor 
with rjv = 1.5 and corresponding value of A'^ on 
the central (lower) full line of FIG. 2 for the upper 
(lower) plot . 



FIG. 4: Normalized VF-spectra (T* G [55,90] 
MeV) for the DE magnetic amplitude. The upper 
and lower figures correspond to the E787 branch- 
ing ratio [13] in TAB. |TI] (the central solid line 
of FIG. 2) and to the NA48/2 branching ratio in 
TAB. HI] (the lower solid line of FIG. 2) respec- 
tively. The solid curves corresponds to a constant 
amplitude, while the dashed and dotted curves 
correspond to form factors with -qv = 0.5 and 
rjv — 1.5 respectively. 



In FIG. 3 we show the normalized r^*-spectra for 
the DE magnetic amplitude with the E787 branch- 
ing ratio [Tg in TAB. [ll] (upper plot) and the 
NA48/2 branching ratio in TAB. [IT] (lower plot). 
The solid curves correspond to a constant ampli- 
tude, the dotted curves to a magnetic form factor 
with rjv = 1-5 and corresponding value of on 



As we can see the changing of the value of the 
branching ratio generates a substantial effect in the 
IF— spectra in FIGS. 5 and 6. 

Subtracting the IB contribution to the W- 
spectrum in ([5]) and assuming E and M constant 



7 



one obtains 



W-spectra : 



dr{E, M) 
dW 



cx INT{E) W'^ + DE{E, M) W'^. (20) 



Then we can fit this to the experimental data de- 
termining E from the interference term, INT(E), 
and M from DE{E,M). Since Ml°^''{z,) in ^ is 
obviously not constant, we would like to question 
whether the presence of this form factor (and no 
interference), could simulate an interference term 
in dini) as observed by NA48/2 [H, [H. In fact 
we plot in FIG. 6 the difference among normalized 



(21) 



1 dr(o,M+(z,)) 



dW 



dr(o,M;°""'-) 

dW 



where the first distribution is generated by the form 
factor structure in eq. (|16|) and the second one by 
constant magnetic term. We study this difference in 
FIG. 6 for DE amplitudes corresponding to form fac- 
tors with rjv = 0.5 (dashed fine) and rjv = 1.5 (dot- 
ted line). Afff and A/'const. are the normalization 
factors for W-spectra, corresponding to the form 
factor and constant magnetic term respectively. 
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FIG. 5: Normalized W-spectra (T; G [0, 80] MeV) 
for the DE magnetic amplitude. The upper and 
lower figures correspond to the E787 branching ra- 
tio [lil in TAB. In] (the central solid line of FIG. 
2) and to the NA48/2 branching ratio in TAB.[ll] 
(the lower solid line of FIG. 2) respectively. The 
solid curves corresponds to a constant amplitude, 
while the dashed and dotted curves correspond to 
form factors with r/v = 0.5 and rjv ~ 1.5 respec- 
tively. 
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FIG. 6: Plots for the quantity A in eq. ((2T|, ob- 
tained by subtracting the normalized W-spectrum 
of constant DE amplitude from the normalized W- 
spectra of DE amplitude corresponding to form 
factors with r/v = 0.5 (dashed line) and r]v = 1-5 
(dotted line); upper and lower plots refer to E787 
branching ratio [l3| in TAB. |lT] (the central solid 
line of FIG. 2) and to the NA48/2 branching ra- 
tio in TAB. |II] respectively (the lower solid line of 
FIG. 2). r<r-range is [0,80] MeV. 



An interference term would appear as a line start- 
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ing from the origin and going to negative values for 
negative interference. It is clear that there is a cor- 
relation between the measurement of the interfer- 
ence and direct emissiom with form factor. We can 
see from FIG. 6 that this depends also on the size 
of the branching ratio: the effect is < 20% for the 
E787 branching ratio and even 30% for the NA48/2 
branching ratio in TAB. |TT] and for the extreme 
value of the VMD parameter, rjv = 1-5 (dotted 
line). However we can exclude that the form factor 
can completely account for the interference effect. 



V. CONCLUSIONS 

VMD has played a very crucial role for our un- 
derstandings of chiral dynamics and its predictive 
power in the strong sector; we expect a similar suc- 
cess in the weak sector, however the path is ap- 
parentely more complicated. The magnetic contri- 
bution to K]^ 7r+7r~7 decays is one of the few 
examples in the weak sector where there is phe- 
nomcnological evidence of VMD. This is particu- 
larly interesting for the interplay with CP violating 
amplitudes and possible New Physics searches. Fur- 
thermore, the precision level of CPT tests is such 
that even A(Kl — > Tr~^T^~j)DE has to be known 
with good accuracy [1^. 

We have studied in this paper the correlated 
channel A'+ — > tt'^tt'^j, in connection with the up- 
coming NA48/2 results. We have analyzed in Sec- 
tion II a set of kinematical variables, Christ's vari- 



ables, used also by NA48/2. 

We have numerically analyzed how to find exper- 
imental evidence of VMD in this channel; compared 
to previous literature [l^ we have payed particular 
attention to the dependence of the parameter ijy 
and , entering in this study, on the measured 
B{K+ ^ TT+n^-f) (see FIG. 2). After this analy- 
sis we have concentrated on the dependence of the 
spectra on VMD see FIGS. 3-5. 

We have explored also the possible dependence of 
the observed interference effect, due to electric tran- 
sitions, on the presence of a form factor in magnetic 
transitions. Since electric transitions are important 
from a chiral dynamics point of view this question 
is relevant: we have found that i) there is a correla- 
tion between the size of the branching ratio and the 
shape of the spectra and ii) for an accurate deter- 
mination of the interference contribution the size of 
the form factor is relevant. 
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